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Transport  Properties  of  Metal  Oxides 
In  the  Presence  of  Sulfur 
H.  S.  Spacil  and  D.  A.  Shores 


I . INTRODUCTION 

The  work  described  in  this  report  has  been  addressed  to  the  effects 
of  sulfur  on  the  oxidation  of  metals  by  gases.  As  distinct  from  the 
type  of  oxidation  known  as  "hot  corrosion",  and  which  involves  formation 
of  condensed  sulfur-containing  salts  on  the  metal  surface,  acceleration 
of  oxidation  rates  of  many  metals  can  occur  in  atmospheres  which  contain 
(>2  but  also  contain  the  oxidized  sulfur-bearing  species  SO2  and  SO^. 

This  acceleration  occurs  when  the  chemical  activity  or  potential  of 
oxygen  in  the  atmosphere  is  sufficient  to  form  the  metal  oxide  MO^  while 
the  activity  or  potential  of  sulfur  is  too  low  to  support  formation  of 
the  sulfide  MS^  from  a thermodynamic  standpoint.*  Two  general  classes 
of  sulfur  effects  have  been  observed  for  metals  which  normally  form 
compact,  protective  oxide  scales  in  C^:  either  a compact  oxide  scale  is 
formed,  but  with  a rate  constant  which  increases  with  sulfur  content  of 
the  atmosphere,  or  a duplex  oxide/  sulfide  scale  is  formed,  often  at  a 
rate  much  higher  than  that  of  simple  oxidation  alone.  In  the  former 


*The  chemical  activities  a^  and  a$  of  oxygen  and  sulfur  are  defined  in 

terms  of  the  equilibrium  partial  pressures,  p_  and  p-  , as  (p.  /p°) 

1 /2  ^ L ^ 
and(p-  /p°)  , respectively,  where  p°  is  the  pressure  chosen  for  the 

2 

reference  state  (usually  101.3  kPa  or  1 atm.).  The  chemical  potentials 
Mq  and  are  defined  as  RT  In  (a^)  and  RT  In  (a$) , respectively,  where 
R is  the  gas  constant  in  Mj/kg-mole/°K  and  T is  the  absolute  temperature 
in  °K.  Units  of  p^  and  p^  are  then  Mj/kg-mole  or  kj/gram-mole. 
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situation,  metal  log raphic  examination  of  the  scale  often  shows  the 
presence  of  sulfide  at  the  metal/scale  interface.  Thus  sulfur  effects 
on  oxidation  can  be  linked  to  an  increase  in  sulfur  activity  or  potential 
to  values  above  those  corresponding  to  thermodynamic  equilibrium  in  the 
atmosphere  to  which  the  metal  is  exposed.  In  general,  oxidation  rates 
of  metals  subject  to  such  effects  are  determined  by  transport  of  metal 
cations  through  the  oxide  scale;  once  a duplex  oxide/sulfide  scale 
forms,  the  much  faster  transport  of  cations  through  the  sulfide  portion 
of  the  scale  would  indeed  be  expected  to  result  in  metal  attack  at  rates 
higher  than  for  simple  oxidation.  The  metals  iron,  nickel,  and  cobalt, 
as  well  as  their  technologically  important  oxidation- resistant  alloys, 
all  are  subject  to  sulfur  effects  in  oxidizing  atmospheres.  Understanding 
the  nature  of  the  increase  in  sulfur  activity  during  oxidation  would 
provide  a basis  for  its  control  or  suppression. 

Prior  work  by  various  investigators  has  suggested  three  possible 
mechanisms  for  observed  sulfur  effects  on  oxidation: 

1)  Surface  reaction  rates  at  the  scale/gas  interface  limit 
the  rate  of  transfer  of  oxygen  from  the  atmosphere  to  the  scale, 
resulting  (under  the  proper  conditions)  in  an  increase  of  sulfur 
activity  at  the  scale  surface;  2)  Gas  phase  transport  through  a 
boundary  layer  of  gas  adjacent  to  the  scale  surface  depletes  the 
oxygen  content  of  gas  within  this  layer,  leading  to  an  increase  of 
sulfur  activity  in  the  gas  at  the  scale  surface;  and  3)  Sulfur 
transport  through  the  scale  itself  either  allows  sulfur  to:  (a)  reach 
a region  within  the  scale  where  the  oxygen  activity  has  decreased 
to  a level  which  allows  sulfide  formation,  or  (b)  results  in  an 





increase  of  sulfur  activity  to  a level  which  likewise  allows  sulfide 
formation. 

Each  of  these  three  mechanisms  will  be  discussed  in  more  detail  in 
Section  IV.  The  present  work  was  concerned  principally  with  the  third 
mechanism.  It  was  originally  proposed  to  measure  the  electrical  con- 
ductivity of  metal  oxides  at  known  oxygen  and  sulfur  activities.  Iron, 
nickel,  cobalt,  and  probably  chromium  oxides  are  metal  deficit,  p-type 
electrical  conductors  in  which  cation  vacancies  are  the  principal 
defect.  The  prototypical  equilibrium  of  cation  vacancies  and  electron 
holes  in  the  oxide  MO  with  oxygen  can  be  written  as 

1/2  o2  <■>:  .r 

where,  in  slightly  modified  Kroger-Vink  notation,  0 L is  a doubly 
(negative)  charged  oxygen  ion  on  an  anion  site  with  zero  excess  charge, 

I 

VM  is  a cation  vacancy  with  a single  negative  excess  charge,  and  h'  is 
an  electron  hole  with  a single  positive  excess  charge.  In  an  intrin- 
sically conducting  oxide,  conductivity  at  a fixed  temperature  is  a 
measure  of  the  highly  mobile  electron  hole  concentration.  If  sulfur 
were  to  enter  the  oxide  with  a non-zero  excess  charge  as  an  electron 
donor  or  acceptor,  a conductivity  change  would  result  from  the  resultant 
change  in  electron  hole  concentration.  Any  such  change  would  also  imply 
a change  in  the  cation  vacancy  concentration,  which  might  explain  the 
increase  of  oxidation  rate  caused  by  sulfur  even  when  no  duplex  oxide/sulfide 
scale  is  formed  during  oxidation. 
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In  practive,  the  absolute  conductivity  need  not  be  measured  to 
determine  relative  conductivity  changes.  The  admittance  Y of  a specimen, 
equal  to  the  ratio  of  current  to  voltage,  can  be  used  to  define  a 
function  <t>'  as 

<f>'  = (y-y  1 ) /y  1 l-l 


where  Y1  is  the  initial  admi  ttance.  Any  change  in  <J>'  is  simply  the 
relative  change  in  conductivity  and  thus  of  electron  hole  concentration. 

Measurements  of  4> ' for  NiO  and  CoO  subjected  to  experimentally 
obtainable  combinations  of  oxygen  and  sulfur  activities  showed  changes 
of  the  order  of  1%  in  the  presence  of  sulfur.  These  measurements,  and 
their  significance,  will  be  discussed  in  detail  in  Section  III.  In  view 
of  the  small  measurable  effect,  the  experimental  program  was  modified  to 
allow  measurement  of  the  defect  diffusivity  of  oxide  specimens  containing 
sulfur.  A step  change  in  0^  partial  pressure  over  an  oxide  causes  a 
corresponding  change  in  the  surface  cation  vacancy  concentration.  As 
vacancies  subsequently  diffuse  in  or  out  of  the  oxide  to  reach  the  new 
equilibrium  value,  both  the  admittance  Y and  the  mass  M of  the  oxide 
will  change  with  time.  The  normalized  change  of  either  is  given  by  the 
function  <J>,  defined  as 


(Y-Yf)/(Y'-Yf) 


(M-Mf)/(M'-Mf) 


1-2 


-5- 


where  the  superscript  i indicates  the  initial  value,  and  as  before,  the 
superscript  f indicates  the  final  value,  and  Y and  M are  the  transient 
admittance  and  mass.  Analysis  of  the  variation  of  <f)  vs.  time  can  yield 
diffusivity  values  for  specimens  of  known  geometry.  These  diffusi vi ties 
are  unambiguous  if  a single  mobile  entity  is  diffusing  via  a single 
process;  changes  in  either  the  diffusion  process  or  the  nature  of  the 
diffusing  entity  lead  to  multiple  diffusivi ties  which  may  be  separable 
if  their  respective  time  constants  differ  sufficiently.  An  advantage  of 
diffusivity  measurements  is  then  the  detection  of  low-level  effects, 
subject  to  the  restriction  of  time  constant  difference.  Such  effects 
were  found  for  both  NiO  and  CoO. 

The  possible  effects  of  sulfur  dissolved  in  a metal  oxide  on  oxida- 
tion rate  and  sulfur  activity  gradients  through  the  oxide  were  also 
considered  theoretically.  Since  electron  holes  and  lattice  defects 
carry  an  excess  charge,  the  local  electrical  potential  must  be  considered 
in  the  formulation  of  appropriate  transport  equations.  This  was  done 
for  a model  of  a growing  oxide  scale,  allowing  two  different  types  of 
sul fur- induced  defects  in  the  oxide.  One  assumed  set  of  defect  proper- 
ties can  be  shown  to  lead  to  increases  in  both  oxidation  rate  and  sulfur 
activity  at  the  metal/scale  interface.  A surprising  implication  of  this 
situation,  however,  is  that  cation  vacancies  diffuse  up  their  own  concen- 
tration gradient,  due  to  the  coupling  between  transport  of  various 
defects  and  electron  holes  through  the  local  electrical  potential. 

Changes  in  cation  vacancy  concentration  of  a metal  oxide  due  to  the 
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presence  of  an  added  component  such  as  sulfur  may  then  not  be  sufficient 
in  themselves  to  explain  changes  in  oxidation  rate. 

I I . Experimental 

A.  NiO 

1)  Polycrystalline  NiO 

Cylindrical  specimens  were  prepared  by  sintering  high  purity  NiO 
powder  in  oxygen  at  temperatures  to  l800°C.  These  had  greater  than  95% 
of  theoretical  density,  and  could  be  ground  to  a diameter  of  0.635  cm 
(0.25  in).  Platinum  electrodes  were  applied  to  allow  4-point  admittance 
measurements. 

A zirconia  electrolyte  EMF  cell  was  utilized  to  monitor  the  instan- 
taneous equilibrium  0^  partial  pressure  over  the  specimen  during  transient 
measurements  in  which  the  admittance  change  with  time  is  obtained 
following  a change  in  nominal  0^  pressure. 

2)  Single  Crystal  NiO 

i A prismatic  single  crystal  about  0.098  cm  on  each  face  by  0.7  cm  in 

length  was  obtained,  with  a total  metallic  impurity  level  of  80  ppm  or 
less*.  The  source  material  had  shown  essentially  intrinsic  conductivity 
behavior  at  high  temperatures,  and  could  thus  be  considered  free  of 
impurities  at  a level  that  would  lead  to  extrinsic  conduction.  Platinum 
electrodes  were  applied  to  allow  4-point  admittance  measurements,  but 
the  EMF  cell  utilized  with  the  polycrystalline  specimens  was  not  employed 
for  this  specimen.  It  was  expected  that  dlffusional  time  constants 


*Atomergic  Chemetals  Corp. , Plainview,  N.J. 
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would  be  much  greater  than  those  for  equilibrium  0 ^ or  partial 

pressure  changes,  and  this  condition  appears  to  have  been  satisfied. 

A series  of  three  runs  at  1000°C  on  this  specimen  will  be  dis- 
cussed. The  description  of  these  runs  is: 

-A 

Run  A.  Specimen  equilibrated  in  CC^  with  0^  pressure  of  10  atm; 

0^  pressure  changed  to  1 atm  and  transient  admittance  measured. 

Run  B.  Specimen  equilibrated  in  02_10%C02;  gas  changed  to  02~ 

10%S02  and  transient  conductivity  measured. 

Run  C.  Specimen  equilibrated  in  02~10$S02;  gas  changed  to  CO2  With 

-4 

02  pressure  of  10  atm  for  105  minutes;  02  pressure  changed  to  1 atm 
and  transient  admittance  measured. 

B.  CoO 

Prismatic  single  crystals  about  0.2  cm  on  each  face  and  from  1.2 

to  1.8  cm  in  length  were  obtained,  with  a total  metallic  impurity  level 

of  380  ppm  or  less*.  The  major  impurity  was  Fe  (200  ppm).  This  impurity 

level  should  ensure  intrinsic  conduction  at  high  temperatures  since  the 

_ “2 

cation  vacancy  concentration  of  CoO  is  of  the  order  of  10  at  1000°C 
for  an  O2  partial  pressure  of  1 atm.  All  measurements  were  made  at 
1000°C  and  1 atm  total  pressure. 

1)  Admittance  Measurements 

Four  different  atmospheres  were  used  to  provide  step  changes  in  O2 

partial  pressures:  02,  90$  O2  - Ar,  air,  and  C02  with  an  O2  partial 

-4 

pressure  of  10  atm.  Transient  admittance  measurements  were  made  for 
five  different  combinations  of  initial  and  final  atmospheres  on  several 
sulfur-free  specimens,  as  shown  in  Table  1 1.1. 

*Atomergic  Chemetals  Corp. , Plainview,  N.J. 
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INITIAL 


90%  02  - Ar 


Table  I I . 1 


FINAL 


90%  02  - Ar  Air  C02 


Transient  admittance  measurements  were  then  made  on  three  crystals  for 
the  conditions  given  in  Table  11.2. 


Table  11.2 


Initial 

Atmosphere 


90%  02  - Ar 
90%  02  - S02 


Equi 1 ibration 
Time 

100  hrs 

50  hrs 


Final 

Atmosphere 
90%  02  - S02 


Specimen 


90%  02  - S02 


50  hrs 


2)  Mass  Measurements 

Three  different  atmospheres  were  used  to  provide  step  changes 
in  02  partial  pressures:  02,  air,  and  Ar  with  an  02  partial  pressure 
estimated  at  10  ^ atm.  Transient  mass  measurements  were  made  for  three 
different  combinations  of  initial  and  final  atmospheres  on  one  specimen, 
as  shown  in  Table  11.3. 
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Table  11.3 


! 


I 


INITIAL  PINAL 

02  Ar 

°2  - X 

Air  X 

Ar  X 

The  same  specimen  was  equilibrated  with  0 2>  then  exposed  to  a 98%  02~ 

S02  atmosphere  for  125  hrs,  and  the  resultant  weight  change  measured. 

A final  transient  mass  measurement  was  made  after  equilibration  with  02 
and  a change  to  an  Ar  atmosphere. 

III.  Experimental  Results 

Mathematical  expressions  for  the  time  dependence  of  the  function, 
ip,  applicable  to  cylindrical  or  prismatic  specimens,  are  well-known^  . 
Where  a single  mobile  entity  is  diffusing  with  a concentrat ion- i ndependent 
diffusivity,  these  expressions  can  be  employed  to  determine  the  diffusivity 
by  numerical  or  graphical  techniques.  For  a metal-deficit  metal  oxide 
in  which  cation  vacancy  diffusion  dominates  the  time  dependence  of  <J>, 
the  measured  diffusivity  will  be  the  chemical  interdiffusion  coefficient 
D.  The  requirement  for  local  electroneutrality  within  the  oxide  leads  to  a 
relation  between  D and  the  vacancy  diffusivity  Dv,  which  can  be  expressed 
as 


Dv  = D/(l  + n) 


1 1 1 - 1 
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where  n is  the  excess  (negative)  charge  on  the  vacancy'  . Further,  the 

self-diffusion  coefficient  for  the  metal,  Du,  is  related  to  Dw  through 

M V 

the  relationship 


j 


(2) 

where  xy  is  the  fraction  of  cation  sites  occupied  by  vacancies 
Independent  measurements  of  D,  DM>  and  n on  several  oxides  have  con- 
firmed these  relations  as  valid  under  conditions  where  the  principal 
defects  are  cation  vacancies.  Since  the  NiO  and  CoO  specimens  used  in 
the  present  work  should  meet  this  criterion,  transient  admittance  or 
mass  measurements  then  yield  D. 

A)  NiO  Admittance  Measurements 
l)  Polycrystalline  NiO 

EMF  cell  voltages  showed  that  the  time  constant  for  an  0^  partial 
pressure  change  in  the  vicinity  of  the  specimen  was  of  the  order  of  100 
seconds,  even  when  changing  from  low  to  high  values.  Such  a change  is 
clearly  not  a step  change,  but  will  not  be  significant  if  the  time 
constant  for  admittance  changes  is  much  greater.  The  time  constant  for 
the  latter  was  of  the  same  order  of  magnitude,  however,  indicating  that 
the  specimen  conductivity  was  following  the  02  partial  pressure  change 
fairly  closely. 
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The  driving  force  for  conductivity  changes  is  the  surface  con- 
centration of  vacancies,  , which  in  turn  will  be  proportional  to  the 

1/(1  + n)  power  of  the  oxygen  activity  at  the  oxide  surface.  A 
normalized  surface  vacancy  concentration  is  given  by  the  function 
$,  defined  as 


I 


$ = 


(s)  ,f 
XV 


111-3 


where  the  superscripts  i and  f have  the  same  meaning  as  in  1-2.  The 
value  of  $ can  be  obtained  readily  from  EMF  cell  voltage  for  any 
appropriate  value  of  n.  For  intrinsically  conducting  NiO,  n will  be 
l for  singly  charged  (negative)  nickel  vacancies  and  2 for  doubly 
charged  nickel  vacancies.  Figure  1 shows  a typical  plot  of  <$>  (solid 
line)  and  $ (dashed  line)  for  polycrystalline  NiO  at  800°C  on  changing 
the  equilibrium  0^  pressure  from  10  to  1 atm.  A value  of  2 was  taken 
for  n to  obtain  but  a value  of  1 does  not  change  the  shape  or  values 
of  $ significantly.  The  points  shown  along  the  solid  line  of  Fig.  1 
are  not  experimental  points,  but  represent  calculated  values  of  <|>  at 
specific  times  as  will  be  discussed  subsequently.  Note  that  the  function 
<t>  lags  the  surface  vacancy  concentration  function  <I>  as  would  be 
expected,  but  that  the  lag  can  be  substantially  less  than  the  total 
time  elapsed.  Thus  a simplified  interpretation  of  normalized  con- 
ductivity changes  to  obtain  D is  not  possible. 


For  a cylindrical  specimen  the  function  <}>  can  be  ex  ressed 
readily  in  terms  of  the  normalized  surface  vacancy  concentration  as 
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oo  f>t 

(J)  = $ - 1|  £ OJ 

k®l 


d$ 

dt  t' 


exp  (-w^  ^-)  dt' 

2 

w, 


I I I -4 


where  t is  the  time  of  measurement,  t'  is  a variable  of  integration, 
and  the  dimensionless  quantity  is  the  root  of  the  zero-order 

Bessel  function  of  the  first  kind.  The  time  constant  x is  the  ratio 

2 ~ (1 ) 

r /D,  where  r is  the  specimen  radius 

Equation  111-4  involves  the  derivative  of  an  experimentally  determined 

quantity,  i.e.,  $.  An  alternative  expression  involving  4>  directly  can 

be  obtained  through  integration  as 


k=l 


, 2 tx 

exp  ' wk  x"' 


t r 


t,  exp  (-wk 


(t-f) 

T 


)df 


m-5 


Since  all  quantities  in  111-5  are  known  except  x,  the  latter  can 
be  treated  as  an  unknown  for  any  pair  of  values  of  <j>  and  t.  This  was 
done  numerically  for  the  data  of  Fig.  1 at  at  0.5  and  t=  134  sec, 
yielding  a value  for  x of  980  sec.  As  a check  on  the  validity  of  this 
value,  it  was  used  to  calculate  <)>  at  the  seven  other  points  shown  on 
Fig.  1.  The  agreement  of  the  calculated  points  with  the  shape  of 
the  experimentally  determined  plot  indicates  that  MI-5  correctly 
describes  the  conductivity  variation  with  time  and  surface  vacancy 
concentration.  The  use  of  an  EMF  cell  can  therefore  allow  determination 
of  time  constants  for  diffusion  even  when  these  are  comparable  to  the 
time  constant  for  the  equilibrium  0 2 pressure  changes  governing  surface 
vacancy  concentrations.  It  may  be  noted  that  the  parameter  n could  also 
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be  treated  as  an  unknown,  thus  in  principle  allowing  a determination  of 

net  vacancy  charge.  In  the  present  treatment  of  the  data,  a value  of  2 

for  n gave  better  agreement  of  calculated  and  measured  conductivities 

than  a value  of  1,  but  the  subject  was  not  investigated  in  detail. 

The  diffusional  time  constant  of  980  sec  and  specimen  radius  of 

“5  2 

0.318  cm  yield  a vacancy  diffusivity  of  3-5  x 10  cm  /sec.  Literature 

“8  2 

data  indicate  that  a value  of  the  order  of  10  cm/sec  might  be  expected 
at  800°C  for  bulk  (single  crystal)  diffusion^.  Either  grain  boundary 
diffusion  or  gas  phase  transport  of  oxygen  through  pores  may  be  responsi- 
ble for  the  high  diffusivity  obtained;  the  latter  is  suggested  as  more 
probable.  The  polycrystalline  NiO  specimens  were  prepared  from  powder 
with  a grain  size  estimated  at  100  pm.  If  diffusion  were  occurring  into 
cylinders  with  a 50  pm  radius,  the  diffusional  time  constant  of  980  sec 
leads  to  a vacancy  diffusivity  of  about  10  an  /sec,  suggesting  that 
this  process  may  actually  be  taking  place.  The  high  apparent  diffusivity 
obtained  appears  in  any  event  to  be  an  artifact  of  the  specimen  structure. 
While  it  may  be  possible  to  correlate  the  apparent  diffusivity  with  bulk 
defect  properties,  this  cannot  be  done  with  certainty.  Attempts  to 
measure  and  characterize  sulfur  effects  on  conductivity  in  polycrystalline 
NiO  are  then  unwarranted  unless  structural  factors  are  better  understood. 

2)  Single  Crystal  NiO 

I 

The  three  runs  made  on  the  single  crystal  specimen  had  these  specific 
objectives.  Run  A was  designed  to  obtain  the  vacancy  diffusivity  in  the 
absence  of  sulfur.  Run  B was  designed  to  maintain  the  0^  pressure  nearly 


i 

* 
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constant, but  to  vary  the  equilibrium  pressure  and  determine  any 
conductivity  changes  due  to  sul fur- induced  defects.  Run  C was  designed 
to  introduce  sulfur  into  the  specimen,  then  obtain  the  vacancy  dif- 
fusivity  while  some  sulfur  still  remained  in  the  specimen. 

Considering  run  B first,  Fig.  2 is  a plot  of  the  relative  con- 
ductivity function  <J>',  as  defined  in  1-1,  vs.  time.  This  figure  shows 
an  initial  decrease  in  conductivity,  followed  by  an  increase  starting 
after  about  25  minutes.  The  conductivity  stabilized  after  several 
hundred  minutes  at  some  value  greater  than  the  initial  value,  then 
underwent  changes  with  relatively  short  time  constants  that  are  most 
probably  associated  with  variations  in  experimental  conditions.  The 
02  - 10%  S02  mixture  was  obtained  by  mixing  the  two  constituent  gases; 
slight  changes  in  the  flow  rate  of  either  would  cause  corresponding 
changes  in  the  02  partial  pressure,  and  thus  in  conductivity.  Flow  rate 
or  compositional  changes  which  affect  gas  mixture  thermal  conductivity  will 
cause  corresponding  changes  in  specimen  temperature  and  in  specimen 
conductivity.  And  the  change  from  an  initial  02  - C02  atmosphere  to  an 
02  - S02  atmosphere  would  also  be  expected  to  change  specimen  temperature 
due  to  the  different  thermal  conductivities  of  these  two  mixtures.  The 
heat  of  activation  for  conduction  in  NiO  is  about  96  kj , so  that  a 1 °C 
temperature  change  of  the  specimen  would  lead  to  relative  conductivity 
change  of  about  0.8%  at  1000°C.  While  Fig.  2 shows  a relative  conductivity 
increase  of  as  much  as  8%,  it  is  estimated  that  variations  in  experimental 
conditions  may  lead  to  relative  conductivity  changes  of  3 to  k%.  Thus 
the  final  relative  conductivity  increase  in  NiO  on  changing  from  02  - 
C02  to  02  - S02  may  be  at  most  of  the  order  of  a few  percent. 
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The  normalized  conductivity  function  of  a prismatic  oxide  specimen 
in  which  cation  vacancy  diffusion  is  occurring  following  a step  change 
in  02  partial  is  given  approximately  by 

In  (<}>)  = In  (a)  - t/x  II  1-6 

in  the  range  of  time  t where  <t>  is  less  than  about  0.6.  In  111-6, 

a is  a numerical  constant  with  a value  of  about  0.8,  and  the  time 

2 2~ 

constant  t is  equal  to  2a  /tt  D,  where  a i*;  one  half  the  prism  face  < 

width,  and  D can  be  calculated. 

Normalized  conductivity  can  be  obtained  for  runs  A and  C in  a 

straightforward  manner  since  conductivity  measurements  could  be  carried 

out  for  sufficient  time  after  an  02  pressure  change  to  obtain  a value 

for  Y^.  For  run  B,  the  value  of  Y^  is  affected  by  the  gas  mixing  system 

as  mentioned  previously.  Thus  the  value  of  Y^  was  taken  as  the  admittance 

reached  during  the  plateau  of  Fig.  2 at  times  of  1*00-500  minutes.  The 

normalized  conductivity  for  run  B will  then  exceed  unity  at  low  times 

during  which  the  conductivity  initially  decreases. 

Figure  3 is  a plot  of  the  normalized  conductivity  function,  on  a 

logarithmic  scale,  against  time  for  the  three  runs.  Run  A shows  the 

expected  linear  (.ortion  with  a time  constant  of  2360  sec  which  corresponds 

“8  2 

to  a vacancy  diffusivity  of  5.6  x 10  cm  /sec,  assuming  a value  for 

n of  2.  While  different  investigators  disagree  on  whether  cation  vacancies 

(3  1*) 

in  NiO  are  singly  or  doubly  charged  , the  mole  fraction  of  these 
is  about  5 x 10  ^ at  1000°C,  leading  to  a value  for  the  self  diffusion 
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of  Ni  as  2.8  x 10  within  the  range  reported^.  The  results  of 
run  A were  quite  reproducible,  and  indicate  diffusion  dominated  by 
cation  vacancies. 

Run  B shows  normalized  conductivity  values  exceeding  unity  at 

low  times,  then  a linear  portion  with  a time  constant  of  3780  sec. 

Although  run  B was  carried  out  at  a nominally  constant  0 2 pressure  while 

S02  was  substituted  for  C02,  the  initial  decrease  in  conductivity 

is  most  probably  due  to  an  02  pressure  change  during  this  substituion. 

The  linear  portion  of  run  B then  represents  diffusion  of  a sulfur- 

induced  defect  which  increases  the  conductivity.  Interpretation  of 

run  B indicates  that  sulfur  entering  NiO  increases  conductivity,  thus 

behaving  as  an  electron  acceptor,  and  induces  a defect  having  an  apparent 
-7  2 

diffusivity  of  10  cm  /sec.  While  interpretation  of  the  linear  portion 
of  run  B is  not  straightforward,  it  is  simple  to  demonstrate  that  the 
normalized  conductivity  function  will  be  dominated  at  long  times  by  the 
slower  of  two  simultaneous  defect  diffusion  processes  provided  the 
conductivity  changes  induced  by  the  different  defects  at  those  times 
are  comparable.  Figure  2 indicates  that  this  latter  criterion  is 
satisfied,  and  thus  that  diffusion  of  the  sul f ur- i nduced  defect  dominates 
run  B after  about  50  minutes. 

Run  C shows  a linear  portion  up  to  about  50  minutes  with  a time 
constant  of  2370  sec  that  can  be  identified  with  cation  vacancy  diffusion. 

A change  in  slope  occurs  after  this  time.  Since  the  specimen  was  held 
at  a low  02  pressure,  after  equilibrating  with  02-10%S02  prior  to  this 
run,  most  of  the  sul fur- induced  defects  would  be  expected  to  diffuse 
out  of  the  specimen  during  this  holding  period  if  the  diffusivity  obtained 
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from  run  B is  correct.  At  long  times,  however,  the  remaining  sulfur- 
induced  defects  might  be  expected  to  affect  the  conductivity.  The 
conductivity  was  increasing  through  run  C,  but  the  removal  of  sulfur- 
induced  defects  which  tend  to  increase  the  conductivity  would  slow  down 
the  rate  of  increase  at  long  times.  Thus  the  direction  of  the  rate  of 
conductivity  change  of  the  second  linear  portion  of  run  C is  consistent 
with  the  interpretation  of  run  B. 

B.  CoO  Admittance  Measurements 

At  an  02  partial  pressure  of  1 atm,  CoO  has  a vacancy  mole  fraction 
of  about  8 x 10  ^ at  1000°C,  so  its  intrinsic  conductivity  should  be 
less  susceptible  to  effects  of  metallic  impurities  than  NiO.  Also,  it 
is  generally  agreed  that,  except  at  higher  temperatures,  singly  charged 
cation  vacancies  are  the  predominant  defect  in  CoO^. 

1)  Oxygen  Activity  Change  in  the  Absence  of  Sulfur 

When  treated  according  to  ill-6,  admittance  data  from  the  runs 
made  as  described  in  Table  I 1.1  of  Section  ll-B-1  yielded  values  of 

D ranging  from  6.38  x 10  ^ to  12.1  x 10  ^ , with  a mean  of  10.0  x 10  ^ 

2 -3 

cm  /sec.  Since  n is  equal  to  1 , a vacancy  mole  fraction  of  8 x 10  3 

leads  to  a value  for  the  self-diffusion  of  Co  at  1 atm  02  partial  pressure 
as  4.0  x 10  This  is  well  within  the  range  reported , and  indi- 
cates diffusion  governed  by  cation  vacancies. 

Results  of  a typical  run  are  shown  in  Fig.  4,  where  the  quantity 
-ln(<)>)  is  plotted  against  time.  The  experimental  points  shown  in 
Fig.  4 give  a satisfactory  fit  to  a straight  line  between  the  limits 

for  <J>  of  about  0.6  and  5 x 10  . The  upper  limit  is  determined  by 
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the  approximat ions  inherent  in  I 1 1-6,  while  the  lower  limit  is  deter- 
minted  by  the  magnitude  of  the  0^  partial  pressure  step  change  together 
with  the  sensitivity  and  stability  of  the  admittance  measuring  cir- 
cuitry. The  lower  limit  corresponds  to  a value  for  -ln(<j>)  of  about 
5;  measurements  out  to  times  leading  to  this  value  are  felt  to  be  reliable. 

2)  Sulfur  Activity  Change  at  Constant  Oxygen  Activity 

Run  8 as  described  in  Table  11.2  of  Section  ll-B-1  allowed 
measurement  of  the  relative  conductivity  change  induced  by  the  intro- 
duction of  sulfur  at  nominally  constant  oxygen  activity.  Results  are 
presented  in  Table  111.1. 


Table  lll.l 


Time 

(seconds) 

CM 

O 

X 

-©- 

0 

0 

300 

0.81 

600 

1.01 

900 

1.42 

1200 

1.52 

1800 

1.63 

21400 

1.73 

3000 

1.94 

3600 

2.04 

4200 

2.25 

4800 

2.25 

5400 

2. 14 

6000 

2.04 

00  - 25000 

2.04  - 2.35 

Table  lll.l  shows  a relative  conductivity  increase  of  about  2$  following 
the  introduction  of  SO^.  The  same  considerations  as  applied  to  run  B 
for  single  crystal  NiO,  as  discussed  in  Section  lll-A-2,  apply  to  run  8 
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for  CoO.  Due  to  a lower  heat  of  activation  for  conduction  in  CoO, 
it  is  estimated  that  variations  in  experimental  conditions  may  lead 
to  relative  conductivity  changes  of  about  1%.  Thus  the  final  relative 
conductivity  increase  in  CoO  on  changing  from  0 2 - Ar  to  02  - S02 
may  be  of  the  order  of  1%. 

3)  Oxygen  Activity  Change  in  the  Presence  of  Sulfur 

Runs  12  and  13  as  described  in  Table  11.2  of  Section  ll-B-1  allowed 
measurement  of  the  normalized  conductivity  function  for  both  positive 
and  negative  step  changes  in  oxygen  activity  following  exposure  of  a 
specimen  to  a sulfur-containing  atmosphere.  Results  of  these  runs  are 
shown  in  Fig.  5;  in  comparison  with  Fig.  4,  it  should  be  noted  that  the 
time  scale  of  Fig.  5 extends  to  higher  values.  An  averaged  curve  is 
shown  for  the  two  runs  in  Fig.  5,  but  each  could  be  separated  approximately 
into  two  near-linear  segments  intersecting  at  a value  for  (f  of  about 
0.025.  This  suggests  that  when  the  initial  conductivity  has  changed 
by  about  98%,  diffusion  of  some  sulfur- induced  defect  determines  further 
conductivity  changes.  Table  III. 2 gives  the  apparent  di ff us i vi t i es 
resulting  from  separation  of  each  curve  into  two  linear  segments:  hese 
are  designated  by  D'  for  short  times  and  D"  for  long  times. 


Table  III. 2 

Run 

pD 1 

2°" 

(cm  ’/sec) 

(cm  /sec) 

12 

9.55  x 10'7 

o 

X 

CM 

13 

10.4  x 10'7 

3.25  x 10"7 
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If  it  is  assumed  that  a sul fur- i nduced  defect  does  increase  the 
conductivity  of  CoO  by  a small  amount  of  the  order  of  1%,  then  Figs.  6a 
and  6b  show  schematically  the  conductivity  variations  with  time  that 
could  give  rise  to  the  results  of  Fig.  5.  In  each  of  these  two  plots, 
the  dashed  lines  labelled  and  are  the  contributions  to  the  total 
conductivity  due  to  electron  holes  associated  with  cation  vacancies  and 
the  sul fur- i nduced  defect,  respectively.  The  solid  line  is  the  total 
conductivity,  which  at  long  times  is  determined  primarily  by  diffusion 
of  the  sul fur- i nduced  defect.  The  apparent  diffusivity  D"  at  long  times 
will  be  less  for  a positive  oxygen  activity  change,  as  in  Fig.  6b  where 
the  individual  conductivity  variations  oppose  each  other,  than  for  a 
negative  oxygen  activity  change,  as  in  Fig.  6a,  where  the  conductivity 
variations  are  in  the  same  direction. 

Run  16  was  essentially  a repeat  of  run  13,  but  without  the  intervening 
exposure  to  02  of  run  12  following  exposure  to  the  90%  02  - S02  atmos- 
phere. The  normalized  conductivity  function  could  again  by  separated 

into  two  near-linear  segments,  yielding  two  apparent  diffusivi ties 

~7  2 

with  D1  having  a value  of  9.28  x 10  ' cm  /sec. 

The  values  of  D1  obtained  in  runs  12,  13  and  16  are  in  good  agree- 
ment with  the  chemical  interdiffusion  coefficient  obtained  in  the 
absence  of  sulfur.  The  short  time  variation  in  the  normalized  con- 
ductivity function  then  can  be  ascribed  to  cation  vacancy  diffusion. 

Since  the  change  to  a lower  apparent  diffusivity  occurs  for  values 
of  -ln(<j>)  less  than  5,  this  change  is  considered  to  be  a reliable 
indication  of  an  increase  of  conductivity  associated  with  a sul fur- i nduced 
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defect  whose  apparent  diffusivity  is  about  5 times  less  than  that  of 
cation  vacancies. 

C.  CoO  Mass  Measurements 

A Cahn  RG  thermobalance  and  digital  data  recording  system  were  used 
to  carry  out  the  measurements  described  in  Section  ll-B-2.  The 
requirement  for  a flowing  atmosphere  over  the  specimen,  however,  caused 
a high  noise  level  in  the  analog  signal  from  the  thermobalance  control 
circuit.  Data  points  taken  by  the  data  recording  system  at  discrete 
times  were  recorded  and  subsequently  computer  processed  by  averaging 
over  an  arbitrary  number  of  adjacent  points  (usually  either  5 or  10). 

This  averaging  removes  some  of  the  information  associated  with  the 
normalized  mass  function  defined  by  1-2,  since  if  carried  to  an  extreme, 
it  will  simply  generate  the  average  value  of  over  the  time  of  the 
experiment.  Values  of  -ln(<}>)  do  not  appear  to  be  reliable  above  3, 
and  results  obtained  are  felt  to  be  less  significant  than  those  obtained 
from  admittance  measurements. 

Step  changes  of  oxygen  activity  in  the  absence  of  sulfur,  as  described 

in  Table  11.3  of  Section  ll-B-2,  resulted  in  averaged  data  yielding 

values  for  the  chemical  interdiffusion  coefficient  of  about  28  x 10  ^ 

2 

cm  /sec  when  treated  according  to  111-6.  The  discrepancy  between  this 

value  and  that  from  admittance  measurements  may  be  due  in  large  part  to 

the  averaging  technique.  Following  exposure  to  a 36%  ' ^0^  atmosphere 

and  then  to  0^,  the  specimen  was  subjected  to  a negative  step  change 

in  oxygen  activity.  An  average  apparent  diffusivity  for  this  run  was 
-7  2 

15  x 10  cm  /sec,  suggesting  that  a slower  diffusion  process  was  influencing 
the  apparent  diffusivity. 
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Thermobalance  data  did  allow  determination  of  the  long  time  mass 

gain  of  the  specimen  after  changing  from  an  atmosphere  to  a 98$ 

C>2  - SO2  atmosphere.  Since  the  rate  of  change  of  mass  is  not  important 

in  such  a measurement,  data  point  averaging  will  improve  the  result. 

Bouyancy  corrections  for  the  difference  in  density  of  the  two  atmospheres 

allow  a determination  of  the  actual  mass  gain  of  the  0.46  gram  specimen 

as  0.094  milligrams.  This  increase  corresponds  to  a sulfur  mole  fraction 
-4 

of  about  5 x 10 

IV.  Discussion 

Other  investigators  have  concluded  that  sulfur  dissolved  in  NiO 
either  decreases  conductivity^  or  has  no  effect  on  conductivity^. 

The  self-diffusion  coefficient  of  sulfur  in  NiO  and  CoO  has  been  reported 
as  about  two  orders  of  magnitude  less  than  the  self  diffusion  coefficient 
for  the  respective  metals^.  The  present  work  indicates  that  exposure 
to  an  atmosphere  containing  10$  S02  increases  the  conductivity  of  NiO 
1 and  CoO  slightly,  at  least  at  high  02  partial  pressures,  and  that 

some  sul fur- i nduced  defect  in  these  oxides  has  an  apparent  diffusivity 
no  more  than  an  order  of  magnitude  less  than  the  chemical  interdiffusion 
coefficient  due  to  cation  vacancies.  For  CoO  at  1000°C  and  1 atm  0,,  the 

; mole  fraction  of  dissolved  sulfur  due  to  exposure  to  an  atmosphere 

r ...  -4 

containing  2$  SO2  is  5 x 10  , or  about  6$  of  the  cation  vacancy  con- 

centration. The  sulfur  activity  of  these  02  - S02  mixtures  is  of  the 
order  of  10  ' ' . Even  at  such  a low  sulfur  activity,  our  results  suggest 
that  a significant  amount  of  sulfur  can  be  dissolved  in  a metal-deficit 
oxide,  but  that  only  a small  fraction  of  this  sulfur  affects  the  electron 
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’ hole  or  vacancy  concentration  at  high  0^  partial  pressures. 

To  consider  whether  the  results  of  the  present  work  are  con- 
sistent with  observed  sulfur  effects  on  oxidation,  it  is  useful  to 
employ  the  reaction  path  concept  of  the  possible  activities  or  chemical 
potentials  that  may  arise  during  oxidation.  This  concept,  as  applied 
to  sulfur  effects  on  oxidation,  utilizes  the  stability  diagram  for 
the  M-O-S  system;  such  a diagram  is  simply  a phase  diagram  depicting 
regions  of  metal,  oxide,  and  sulfide  stability  in  terms  of  oxygen  and 
sulfur  activities  or  chemical  potential  rather  than  composition. 

The  locus  of  all  possible  simultaneous  values  of  oxygen  and  sulfur 
activities  encountered  during  oxidation  (the  "reaction  path"  is  super- 
imposed on  the  stability  diagram  to  explain  or  predict  the  existence 
of  observed  phases.  Reaction  paths  have  been  presented  by  several 
, investigators  in  connection  with  all  three  of  the  mechanisms  for 

sulfur  effects  on  oxidation  summarized  in  Section  | (9.10,11)  |n 
general,  as  stated  in  Section  I,  the  objective  of  such  an  approach  is 
to  ascertain  whether  increases  in  sulfur  activity  can  occur  during 
metal  oxidation. 

A.  Surface  Reaction  Rate 

The  reaction  path  concept  is  least  helpful  here  since  surface 

reaction  rate  limitation  of  scale-forming  reactions  implies  discontinuities 

in  oxygen  or  sulfur  activities,  or  both,  at  the  scale/gas  interface. 

Surface  reaction  rate  limitation  of  sulfide  formation  on  Fe  has  been 

demonstrated,  but  only  in  atmospheres  where  sulfides  were  thermo- 

(12) 

dynamically  stable  . A slow  adsorption  process  at  the  scale  surface, 

. 
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for  example,  can  lead  to  surface  reaction  rate  limitation  of  scale 
formation.  For  formation  of  duplex  oxide/sulfide  scale,  surface  reaction 
rates  would  have  to  affect  both  oxygen  and  sulfide  activities  at  the 
scale  surface,  especially  in  atmospheres  with  a high  0 2 partial 
pressure:  the  oxygen  activity  would  have  to  decrease,  while  the  sulfur 

activity  would  have  to  simultaneously  increase.  Since  duplex  scale 
formation  is  observed  on  several  different  metals  and  alloys,  this 
characteristic  of  surface  reaction  rates  would  have  to  apply  to  a variety 
of  oxide  or  oxide/sulfide  surfaces.  In  view  of  the  highly  specific 
nature  of  surface  reaction  rate  kinetics,  this  wide  applicability  seems 
unlikely. 

B.  Gas  Phase  Transport 

Any  solid  metal  exposed  to  a flowing  atmosphere  will  have  a mass 
transport  boundary  layer  of  gas  adjacent  to  its  surface.  The  gas 
velocity  will  decrease  throughout  this  layer,  becoming  essentially 
zero  at  the  solid  surface.  Diffusive  transport  of  oxygen  and  sulfur 

< 

must  occur  through  this  layer  for  oxide  and  sulfide  to  form.  If  oxygen 
is  removed  from  the  gas  at  a sufficiently  high  rate,  composition 
changes  in  the  boundary  layer  can  cause  the  sulfur  activity  of  the  gas 
at  the  metal  surface  to  increase  while  the  oxygen  activity  simultaneously 

J decreases.  Figure  7 shows  the  reaction  path  for  such  a situation;  a 

metal -oxygen-sul fur  system  in  which  only  metal,  oxide  MO,  and  sulfide 
MS  can  exist  has  been  assumed  for  simplicity.  The  coordinates  are  the 
chemical  potentials  pg  and  p^  of  oxygen  and  sulfur.  The  stability  fields 
for  each  possible  phase  are  indicated,  and  are  separated  by  straight 
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solid  lines  which  represent  the  chemical  potentials  at  which  two  phases 
can  co-exist.  The  bulk  gas  composition  has  been  taken  to  lie  in  the  MO 
stability  field.  Diffusive  transport  of  oxygen  and  sulfur  through  the 
boundary  layer  may,  under  the  appropriate  conditions,  produce  the 
reaction  path  shown,  leading  to  thermodynamic  stability  of  both  MO  and 
MS  at  the  scale  surface.  Duplex  scale  formation  would  be  expected  under 
such  conditions. 

Many  laboratory  investigations  of  oxidation  by  oxygen  and  sulfur 
have  been  carried  out  by  using  gas  mixtures  of  SC^  and  some  inert  gas 
such  as  Ar.  Duplex  oxide/sulfide  scale  formation  can  be  achieved  readily 
over  a range  of  SO2  concentrations  and  temperatures,  but  tends  to  reach 
a maximum  rate  for  an  intermediate  SO^  concentration.  To  determine  if 
this  result  were  likely  to  be  a consequence  of  known  factors  rather  than 
some  property  of  the  oxide  or  sulfide,  the  diffusive  transport  of  oxygen 
and  sulfur  through  the  boundary  layer  was  examined.  In  particular,  the 
conditions  for  a sulfur  activity  increase  leading  to  sulfide  formation 
at  the  surface  of  an  existing  oxide  scale  were  considered  since  a pre- 
existing oxide  layer  will  prevent  duplex  scale  formation  if  sufficiently 
thick. 

Figure  8 shows  schematically  the  gas  phase  compositions  through 
the  boundary  layer  of  thickness  6 adjacent  to  metal  covered  with  an 
oxide  layer  of  thickness  L for  the  condition  where  net  sulfur  transport 
is  zero,  but  the  sulfur  activity  has  just  risen  to  the  point  where  sulfide 
will  be  stable  at  the  oxide  surface.  The  SO2  and  contents  of  the 
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bulk  gas  are  the  mole  fractions  xb_  and  x|?  , respectively.  This  0_ 

2 z 1 

content  of  the  bulk  gas  in  a laboratory  experiment  will  depend  on  the 

purities  of  the  S02  and  inert  gas  used;  while  the  inert  gas  can  be 

purified  to  very  low  0^  levels,  the  S02  will  contain  some  residual  0 2 in 

( 1 3) 

amounts  up  to  1000  ppm  . At  the  oxide/gas  interface,  the  S02 
content  of  the  gas  is  essentially  unchanged,  while  the  02  content  has 
dropped  to  Xq  and  the  S2  content  has  rise  from  a near-zero  value  in  the 
bulk  gas  to  a value  of  x!  . Cation  transport  of  metal  through  the  oxide 
to  the  oxide/gas  interface  where  metal-oxygen  reaction  actually  occurs 
is  responsible  for  the  gradient  in  02  content. 

In  order  for  sulfide  to  be  stable  at  the  oxide  surface,  the  oxide 
thickness  must  be  less  than  some  critical  value,  denoted  by  Lc<  Thicker 
oxide  layers  result  in  cation  transport  which  is  too  low  to  support  the 
required  decrease  in  02  content  that  will  allow  the  sulfur  activity  to 
increase  at  the  oxide  surface.  For  cation  transport  via  singly  charged 
vacancies,  a straightforward  mass  transport  balance,  together  with  the 
thermodynamic  condition  for  simultaneous  oxide/sulfide  co-existence, 
leads  to  an  expression  for  the  ratio  l_c/5  as 


Lc/6~ 


J/6,  b x 1/6  , m1/2 

K (xso2}  " (a0} 

xb  - K2/3(xb  )2/3 

o2  K uso2; 


where  the  constant  K is  given  by  K^S/KM0KS0  ’ tbe  equ‘ 1 ibrium  constants 
for  formation  of  MS,  MO,  and  S02  from  the  elements,  and  ag  is  the  oxygen 
activity  at  the  metal/oxide  interface,  or  1/K^q-  The  proportionality 
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constant  implied  by  IV-1  contains  the  appropriate  d i f f us i v i t i as  for  the 
oxide  and  gas  phases.  It  is  the  variation  of  the  ratio  L c /6  with  bulk 
gas  S02  content  that  is  of  interest,  however. 

Figure  9 is  a normalized  plot  of  L^/6  obtained  by  assuming  that  the 
inert  gas  and  the  S02  contain  1 and  1000  ppm  of  02,  respectively,  and  by 
utilizing  the  equilibrium  constants  for  formation  of  CoS,  CoO,  and  S02 
at  800°C.  This  figure  indicates  that  at  both  very  low  and  very  high  S02 
contents  of  the  gas,  a relatively  thin  layer  of  pre-existing  oxide  will 
prevent  sulfide  formation  at  the  oxide  surface.  Intermediate  S02 
contents,  however,  require  the  thickest  preexisting  oxide  layer  to 
suppress  formation  of  duplex  oxide/sulfide  scales.  Once  sulfide  sta- 
bility at  the  oxide  surface  has  been  reached,  the  rates  of  formation  of 
oxide  and  sulfide,  and  thus  the  amounts  of  each  formed,  will  be  determined 
by  metal  transport  in  the  sulfide  as  well  as  the  oxide,  and  the  preceding 
analysis  no  longer  applies.  It  is  thus  concluded  that  initiation  of 
duplex  oxide/sulfide  scale  formation  can  be  accounted  for  by  expected 
transport  and  thermodynamic  porperties  of  the  oxide  and  gas  phases. 

C.  Sulfur  Transport  Through  Scale 

While  gas  phase  transport  conditions  may  lead  to  duplex  oxide/sulfide 
> scale  formation,  an  increase  in  sulfur  activity  through  an  oxide  scale 

is  required  to  explain  formation  of  sulfide  at  a metal/oxide  interface  in 
an  atmosphere  which  has  a sulfur  activity  below  that  of  the  metal/sulfide 
equilibrium.  This  situation  would  lead  to  a reaction  path  such  as  shown 
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in  Fig.  iO,  where  a negligible  change  in  oxygen  and  sulfur  activities 

across  the  gas  phase  boundary  layer  has  been  assumed.  Reaction  paths 

for  activities  across  a scale  in  the  simultaneous  presence  of  oxygen  and 

sulfur  have  usually  been  presented  with  the  sulfur  activity  or  potential 

decreasing  toward  the  metal/scale  interface;  an  exception  is  for  a scale 

which  may  contain  pores^'^.  The  gas  reaching  such  pores,  through 

fissures  in  the  scale  which  are  presumably  associated  with  the  pores, 

loses  oxygen  with  an  accompanying  rise  in  sulfur  activity.  An  analogous 

mechanism  for  a rise  in  carbon  activity  across  the  oxide  scales  formed 

(14) 

during  oxidation  of  Fe  by  CC^  has  been  demonstrated  . There  is 
little  evidence  that  oxide  scales  formed  in  0 2 - SO^  atmospheres  contain 
pores  which  could  similarly  cause  an  increase  in  sulfur  activity  across 
the  scale.  And  if  only  narrow  fissures  through  the  scale  are  present, 
removal  of  oxygen  from  the  gas  penetrating  these  implies  their  closure 
by  oxide  formation.  Further,  an  increase  in  sulfur  activity  across  a 
scale  due  to  gas  reaching  pores  does  not  account  for  increases  observed 
in  the  oxidation  rate  even  when  no  duplex  oxide/sulfide  scale  is  present^"^ 
It  is  then  appropriate  to  consider  the  related  questions  of  whether 
sul fur- induced  defects  in  a growing  oxide  can  lead  to  an  increase  in 
sulfur  activity  at  the  metal/oxide  interface,  and  also  lead  to  an  increase 
in  oxidation  rate  without  the  formation  of  duplex  scales.  Rather  than 
select  a particular  oxide,  the  prototypical  metal  deficit  oxide  MO 
having  properties  representative  of  real  oxides  will  be  examined. 
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1)  Oxidation  Rate  in  the  Absence  of  Sulfur 
It  will  be  assumed  that  the  principal  defects  in  the  oxide  MO 
are  singly  charged  metal  vacancies.  The  concentration  of  these 
at  an  oxide  surface  exposed  to  some  fixed  oxygen  activity  is  then 
determined  by  the  equilibrium 

1/2  02(g)  t 0=  | o + Vm  + h’  ,V"a 

where  the  symbols  have  the  same  significance  as  for  the  prototypical 
equilibrium  discussed  in  Section  I.  Letting  denote  the  equilibrium 
constant  for  IV-a,  the  concentrations  of  cation  vacancies  and  electron 
holes  are  related  by 


K0  = 


Vh/ao 


I V-2 


where  x„  and  x,  are  the  mole  fractions  of  vacancies  and  electron  holes, 
V h 

respectively. 

Only  cation  vacancies  and  electron  holes  are  considered  to  have 
an  appreciable  mobility  in  the  oxide.  The  transport  equations  for 
the  fluxes  of  electron  holes  and  vacancies,  and  Tv,  respectively, 
through  an  oxide  of  thickness  L can  be  written  as 


L 


c 


dlnx. 

■<Vh  <RT  — 


iv-3 
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and 


dlnXy 

Vv  (RT  “dX 


C dUx 

F dP 


IV-** 


where  X is  a normalized  distance  through  the  oxide  (distance  divided 
by  L) , 8^  and  8y  are  the  mobilities  of  electron  holes  and  vacancies, 
respectively,  U is  the  local  electrical  potential,  T is  the  absolute 
temperature,  R is  the  gas  constant,  F is  the  Faraday,  and  c is  the 
concentration  of  cation  sites  in  the  oxide.  With  this  formulation, 
the  mobility  of  any  entity  is  proportional  to  its  diffusivity;  the 
Nernst-Einstein  relation  gives  the  proportionality  constant  as  1/RT. 

Two  conditions  applicable  to  the  oxide  are  the  requirements  for 
local  electroneutrality  and  for  zero  net  current  density.  These  con 
ditions  can  be  expressed  as 


- x„ 


I V- 5 


and 


respectively,  for  the  prototypical  oxide. 

The  three  equations  I V- 2 , I V- 5 and  IV-6  can  be  used  to  obtain 

three  simultaneous  first  order  ordinary  differential  equations  in 

the  variables  x.  , xw,  and  U.  While  I V- 2 and  I V- 5 can  be  differentiated 
h V 

directly,  I V- 3 and  I V- 4 can  be  substituted  into  IV-6  to  obtain  the 


i 
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differential  form.  The  three  differential  equations  are 

dlnx,  + dlnx.,  = dlna. 
h v u 

x.dlnx,  - x dlnx..  = 0 

h h v V 

-RT  B,  x.dlnx,  + RT  B,.x  dlnx.  - F(B.x  + Bwx.,)dll  = 0 
hh  n vv  V hh  VV 

Equations  IV~7a  through  IV-7c  allow  evaluation  of  the  derivatives 

dlnx./dlna.,  dlnx./dlna..,  and  dU/dlnaft  at  any  values  of  x.  and  x... 
n u v u u n v 

Since  U does  not  appear  explicitly  in  these  equations,  the  absolute 
value  of  the  electrical  potential  cannot  be  determined,  and  U may  be 
given  any  convenient  value  at  some  corresponding  value  of  aQ. 

For  a growing  oxide,  a^  will  be  used  to  denote  the  oxygen  activity 
in  the  bulk  gas  to  which  the  oxide  is  exposed,  and  a*  is  the  oxygen 
activity  at  the  metal/oxide  interface.  Calculation  of  the  values  of  x^, 
Xy,  and  U at  any  point  within  the  oxide  then  becomes  an  initial  value 
integration  problem:  I V- 2 and  I V- 5 are  used  to  obtain  the  values  of  xh 
and  Xy  at  the  oxide/gas  interface  where  the  oxygen  activity  is  a^,  U 
can  be  assigned  a value  of  zero  at  this  interface,  and  integration  of 
the  three  above  derivatives  is  carried  out  from  a^  to  any  desired  oxygen 

it 

activity  greater  than  or  equal  to  a^.  It  is  convenient  to  let  the 
quantity  RT/F,  with  the  dimensions  of  volts,  be  represented  by  ip;  then 
U can  be  expressed  as  the  dimensionless  ratio  U/ip. 

The  metal  flux  through  the  growing  oxide  is  equal  in  magnitude  to 
the  vacancy  flux,  but  with  the  opposite  sign.  Since  there  are  no 
sources  or  sinks  of  metal  within  the  oxide,  the  metal  flux  must  be 
constant  across  the  oxide,  and  independent  of  the  normalized  distance  X. 


IV- 7a 
IV-7b 
IV-7c 
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A quantity  D,  with  the  dimensions  of  a diffusivity,  can  then  be  defined 
as 


D = - 


IV- 8a 


or,  from  IV-A,  and  assuming  that  8y  is  constant  across  the  oxide, 

b 

P dmx  . ... 

3T^)d,nao  ,''-8b 

a0 

It  is  apparent  that  D,  which  is  directly  proportional  to  the  metal 
flux  and  inversely  proportional  to  the  oxide  thickness,  wiii  be  related 
to  a rate  constant  for  oxidation.  Since  D is  independent  of  X,  it 
must  also  be  independent  of  a^.  Then  in  principle  its  value  could  be 
obtained  at  any  single  point  within  the  oxide.  The  integral  formulation 
of  I V- 8 avoids  the  necessity  for  selecting  such  a point. 

This  prototypical  example  actually  can  be  solved  in  closed  form, 
and  leads  to  the  well-known  situation  of  parabolic  oxidation.  The 
quantity  D is  given  by 


D - l/2(WN0/PM0W0)2kp 


I V-8c 


where  and  are  the  molecular  weights  of  MO  and  0,  respectively, 

Pun  is  the  density  of  MO,  and  k is  the  parabolic  rate  constant  for 
MU  p 

weight  gain.  0 can  also  be  expressed  as 


2BTBvKi/Z 


[<ab),/2 


<.;>,/2i 
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for  the  condition  that  » 6^.  More  complex  situations,  however, 
can  be  treated  by  the  same  approach  of  formulating  the  appropriate 
differential  equations  and  integrating  over  a^.  If  metal  transport 
takes  place  only  by  vacancy  motion,  equation  I V- 8b  will  still  yield 
the  correct  value  for  D.  In  particular,  if  additional  defects  which 
do  not  cause  metal  transport  are  to  be  considered,  then  the  value  of 
D can  be  compared  to  a reference  value  D°  evaluated  in  the  absence 
of  such  defects  and  integrated  over  the  limits  of  0 and  1.  Thus  for  the 
prototypical  oxide  MO,  D°  is  given  by 

D°  = 2RT  S^2  IV-9 

again  subject  to  the  condition  that  6^  » 

The  variation  of  x,  , x„,  and  U across  the  oxide  MO  are  shown  in 
n v 

“8  b 

Fig.  11,  where  Kq  has  been  taken  as  10  and  a“  is  unity  (1  atm  02) . 

The  cation  vacancy  and  electron  hole  concentrations  are  equal,  to  satisfy 
the  electroneutrality  and  zero  net  current  conditions.  It  was  assumed 
that  electron  hole  mobility  was  much  greater  than  vacancy  mobility. 

While  both  electron  holes  and  vacancies  flow  down  their  respective 
concentration  gradients,  the  electrical  potential  assists  the  motion 
of  negatively  charged  vacancies  and  opposes  the  motion  of  positive 
electron  holes. 

2)  Oxidation  Rate  in  the  Presence  of  Sulfur 

It  might  be  expected  that  the  electrical  potential  would  affect 
the  distribution  of  sul fur- i nduced  defects.  The  present  work  indicates 
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that  at  high  partial  pressures,  such  defects  do  not  greatly  influence 
oxide  conductivity.  Thus  the  majority  of  these  cannot  carry  any  excess 
charge.  Other  work,  however,  suggests  that  sulfur  enters  metal  oxides 
as  an  electron  donor  to  reduce  the  electron  hole  concentration;  this 
behavior  may  reflect  a sul fur- induced  defect  present  at  low  02  partial 
pressures ^ . Therefore  two  types  of  su  1 fur- i nduced  defect  in 
oxide  MO  were  postulated. 

The  first  defect  is  the  type  suggested  by  this  study,  and  consists 
of  S++  ions  on  cation  sites.  The  equilibrium  governing  these  is 

1/2  S2(g)  + V'  + h-  l IV-b 

This  defect,  with  a concentration  xj,  would  have  no  effect  on  electrical 
conductivity  (since  the  hole  concentration  is  still  determined  by  I V- 2 ) , 
and  is  assumed  to  have  a constant  mobility  of  10  times  less  than  the 
vacancy  mobility.  The  ionic  radius  of  S ions,  which  would  determine 
whether  such  defects  are  physically  reasonable,  is  difficult  to  evaluate 
since  few  solid  compounds  contain  this  ion.  The  ionic  radius  of  sulfur 
decreases  in  going  from  S°  to  S , and  it  would  be  expected  that  S 
follows  this  trend.  While  the  experimental  observations  of  the  present 
work  regarding  mobility  of  sul fur- induced  defects  at  high  0^  pressure 
depends  upon  a slight  increase  in  conductivity  due  to  these  defects,  it 
also  appears  that  only  a small  fraction  of  such  defects  are  affecting  the 
conductivity.  Thus  the  postulated  defect  has  zero  excess  charge. 


-35- 


The  second  defect  has  been  suggested  previously  as  noted,  and 
consists  of  S ions  on  anion  sites.  The  equilibrium  governing  these  is 

1/2  S2(g)  + 0=|  q + h-  ^ S"  | q + 1/2  02(g)  IV-c 

This  defect,  with  a concentration  x'^,  would  decrease  electrical  conductivity, 

especially  at  low  0^  partial  pressures,  and  is  assumed  to  have  a constant 

i, 

mobility  10  times  less  than  the  vacancy  mobility.  This  mobility  is 
based  on  the  low  self-diffusion  coefficient  observed  for  oxygen  in 
metal-deficit  oxides. 

The  equilibrium  constants  for  the  two  postulated  defects  are 

Ks  ■ Xs  'V/h  lv'10 

. j-  ^ Ks  - x'sWh  lv-" 

i Equations  IV-10  and  I V- 1 1 can  be  differentiated  directly.  The  fluxes 

1 1 

and  for  the  two  sul  fur- induced  defects  can  be  expressed  in  a form 
similar  to  that  of  \\l-k.  The  conditions  of  local  electroneutrality  and 

zero  net  current  density  still  apply  to  the  oxide  containing  sulfur- 

i 

j induced  defects.  An  independent  specification  of  the  sulfur  flux  is 

> 

required  for  complete  definition  of  the  problem,  however.  An  appropriate 
condition  is  that  the  total  sulfur  flux,  the  sum  of  and  r^',  be  zero. 

If  the  sulfur  activity  at  the  metal/oxide  interface  is  less  than  the 
value  required  to  form  sulfide,  then  the  metal  will  not  constitute  a 
sink  for  sulfur  and  no  sulfur  transport  will  occur.  Imposing  a zero 
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net  sulfur  flux  will  allow  determination  of  the  sulfur  activity  that  can 
be  attained.  Once  sulfide  stability  is  reached,  a net  sulfur  flow  will 
occur  and  a different  anlysis  would  have  to  be  considered  to  determine 
the  rate  of  sulfide  formation.  But  with  the  zero  net  sulfur  flux  specified, 
six  simultaneous  first  order  ordinary  differential  equations  analogous 
to  I V- 7 can  be  obtained  in  the  variables  xh,  xy,  U,  x^ , x1^,  and  a$. 
Calculation  of  these  again  becomes  an  initial  value  integration  problem, 
with  the  integration  taking  place  over  oxygen  activity  from  a^  to 
a^.  Neglecting  any  possible  metal  transport  by  the  first  type  of  sulfur- 
induced  defect,  the  value  of  D can  still  be  obtained  from  I V- 8b . The 
differential  equations  can  no  longer  be  solved  in  closed  form,  and 
a numerical  integration  procedure  capable  of  handling  moderately  stiff 
equations  must  be  employed. 

Table  IV. 1 presents  the  results  of  solutions  obtained  at  several 
values  of  the  bulk  sulfur  activity  a|?  (which  is  an  initial  condition 
for  the  integration).  Other  initial  conditions  at  the  oxide/gas  inter- 
face are:  a^  = 0. 1 and  values  of  x^  = 1.3  x 10  ^ and  xj.'  = 2.8  x 10  ^ 
for  the  specific  case  of  a^  = 10  These  latter  two  conditions  are 

consistent  with  the  assumed  properties  of  the  sul fur- induced  defects, 
while  the  bulk  gas  oxygen  and  sulfur  activities  selected  for  calculations 
represent  a gas  with  1%  02  and  an  S02  content  ranging  from  very  low  to 
about  10%.  In  Table  IV. 1,  x^,  x^,  and  L)1  are  the  values  of  these  quantities 

at  the  metal/oxide  interface,  and  D/D°  is  the  ratio  of  D to  D°  calculated 

-8 

from  I V- 9 • The  value  of  Kg  was  again  taken  as  10  , leading  to  a 

vacancy  concentrat ion  at  the  oxide/gas  interface  of  3.16  x 10  ■*. 
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Table  IV. 1 


b 

as 

1* 

aS 

X * 

V 

u'/ty 

D/O0 

IQ'18 

2. 0A 

X 

io'19 

1.00  X 

IO'8 

7-90 

0.319 

io-'6 

2.0U 

X 

io*17 

8.91  x 

IO*9 

7.79 

0.319 

io-^ 

2.51 

X 

io'15 

1.29  x 

io'9 

5.87 

0.322 

IO'13 

8.71 

X 

io'1 1 

2.88  x 

io-6 

6.85 

0.31*7 

IO-12 

3.63 

X 

O 

l 

oo 

5-75  x 

10~5 

7.61 

0.370 

io'n 

3-89 

X 

10'6 

5.89  x 

io-4 

8.72 

0.1*1 1 

io’10 

2.1*5 

X 

O 

I 

XT 

A. 68  x 

io'3 

9.63 

0.1*87 

One  feature  of  Table  IV. 1 is  immediately  apparent:  the  postulated 
defects  lead  to  an  increase  in  metal/oxide  interface  sulfur  activity 

- 1 l* 

over  that  of  the  bulk  gas  at  sulfur  activities  exceeding  about  5 x 10  , 

, as  well  as  to  an  increase  in  parabolic  oxidation  rate  constant  (as 

reflected  by  D/D° , recognizing  that  D/D°  = 0.319  for  the  limiting  case 

of  a^  = 0,  a^  = 0.1).  It  is  also  apparent  that  these  results  are  not 

simply  explicable.  The  electrical  potential  at  the  metal/oxide  interface 

first  decreases  with  bulk  sulfur  activities,  and  then  increases.  At 

- 1 2 

bulk  sulfur  activities  of  10  or  higher,  the  vacancy  concentration 
at  the  metal/oxide  interface  exceeds  that  at  the  oxide/gas  interface. 

Thus  cation  vacancies  are  diffusing  up  their  concentration  gradient, 
with  a resulting  parabolic  oxidation  rate  constant  that  is  higher  than 
if  they  had  been  diffusing  in  a "normal"  manner.  Presumably  the 
electrical  potential  still  affects  vacancy  motion,  but  a conclusion 
is  that  multiple  defect  interactions  must  be  treated  as  rigorously  as 

i # 

!♦ 
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possible  to  obtain  correct  results  regarding  effects  on  chemistry  and 
kinetics. 

The  conditions  leading  to  the  results  of  Table  IV. 1 were  not 
chosen  to  obtain  predetermined  values.  The  sulfur-induced  defects  were 
modelled  in  view  of  existing  experimental  evidence.  Different  concentrations 
and  mobilities  of  the  su 1 fur- i nduced  defects  lead  to  somewhat  different 
but  generally  similar  results.  Other  types  of  sul fur- i nduced  defects 
can  be  postulated,  of  course,  leading  to  completely  different  results. 

The  success  of  the  postulated  sul fur- i nduced  defects  at  accounting  for 
an  increase  in  sulfur  activity  across  a growing  oxide  and  an  increase 
in  parabolic  oxidation  cannot  be  taken  to  mean  that  such  defects  exist. 

The  prototypical  oxide  examined  here  has  in  effect  been  considered 
as  a single  crystal;  grain  boundary  defects  may  well  overwhelm  any 
lattice  defects,  though  the  former  could  presumably  lead  to  many  of  the 
same  results  as  the  postulated  lattice  defects.  This  examination  of 
possible  defect  interactions  in  a growing  oxide  exposed  to  a sulfur- 
containing  atmosphere  does  show  that  values  of  sulfur  activity  or 
potential  along  the  reaction  path  across  the  oxide  cannot  be  assumed 
to  simply  remain  equal  to  the  bulk  gas  value  or  even  decrease. 

V.  Conclusions 

1)  At  high  0^  partial  pressures,  exposure  of  intr inscial ly 
conducting  single  crystal  NiO  and  CoO  to  SO 2 results  in  small  conduc- 
tivity increases  of  the  order  of  1$  at  1000°C. 

2)  The  apparent  diffusivity  of  defects  in  single  crystal  NiO 
and  CoO  which  have  been  exposed  to  S02,  as  measured  by  conductivity 
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* variation  with  time  following  a change  in  0^  partial  pressure,  shows 

a decrease  at  long  times  that  appears  to  be  associated  with  a 
sulfur-induced  defect. 

I 

3)  Thermogravimetric  measurements  indicate  a sulfur  solubility 
of  about  0.05  mole  percent  in  CoO  exposed  to  0 2 " 2$  S02  at  1000°C. 

b)  Theroetical  considerations  show  that  the  presence  of  a 
lattice  defect  induced  by  sulfur,  and  having  properties  consistent  with 
the  above  observations,  can  lead  to  an  increase  in  sulfur  activity 
or  chemical  potential  across  a growing  layer  of  metal-deficit  oxide. 

The  above  conslusions  suggest  that,  at  high  02  partial  pressures, 
the  principal  sulfur- induced  defect  in  divalent  metal-deficit  oxide 
, may  be  an  S ion  occupying  a cation  site.  The  existence  of  such  a 

defect  would  net  rule  out  the  existence  of  other  sul fur- induced  defects 

I 

at  low  02  partial  pressures. 

i 

I 
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Fig.  1 Normalized  conductivity  (<p,  solid  line)  and  normalized  surface  vacancy  concen 
(t,  dashed  line)  for  polycrystal  1 ine  NiO  at  800  °C.  Circles  are  calculated  values. 
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Fig.  3 Normalized  conductivity  for  single  crystal  NiO  at  1000  °C. 
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Fig.  6 Schematic  variation  of  conductivity  with 
time  after  step  change  in  O2  partial  pressure 
for  oxide  crystal  containing  sulfur-induced 
defects . 


Fig.  7 Reaction  path  for  formation  of  duplex  oxide/sulfide  scales  by  in 
crease  of  sulfur  potential  in  gas  phase  boundary  layer. 
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thickness  for  formation  of  duplex  oxide/sulfide  scales 


